Enhanced Superconductivity in Superlattices of High-T c Cuprates 
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The electronic properties of multilayers of strongly correlated models for cuprate superconduc- 
tors are investigated using cluster dynamical mean- field techniques. We focus on combinations of 
underdoped and overdoped layers and find that the superconducting order parameter in the over- 
doped layers is enhanced by the proximity effect of the strong pairing scale originating from the 
underdoped layers. The enhanced order parameter can even exceed the maximum value in uniform 
systems. This behavior is well reproduced in slave-boson mean-field calculations which also find 
higher transition temperatures than in the uniform system. 
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The fabrication and characterization of artificial het- 
erostructures involving metallic-ion oxides is one of the 
main topics of current materials science 0, 0, 01 ■ Het- 
erostructure might be useful to explore new and bet- 
ter functionalities that cannot be achieved in the bulk 
form. Among various oxides.of particular interest are 
correlated-electron materials [4] which are known to ex- 
hibit rich phenomena including high-T c superconductiv- 
ity in the cuprates [5fl . 

In this Letter, we are interested in the properties of 
model superlattices (SL's) for high-T c cuprates with dif- 
ferent doping levels, in particular, a combination of un- 
derdoped (UD) and overdoped (OD) cuprates. In the for- 
mer, a high pairing scale is observed both experimentally 
@ and theoretically 0, while in the latter, stronger co- 
herence of carriers exists, although T c is lower than at op- 
timal doping in both the regions. By connecting the two 
doping regions in a multilayer geometry, the proximity of 
the strong pairing interaction from the UD region may 
produce stronger superconductivity in the OD region. 
Similar proposals are seen in Refs. d, Conversely, 
the phase stiffness of a superconductor with small super- 
fluid density may be increased by coupling it to a good 
metal resulting in a higher T c [lfjj. The effects of spatial 
inhomogeneity in (quasi-) two-dimensional systems have 
been investigated in the light of the nanoscale inhomo- 
geneity often observed in UD cuprates 11, 12, 0, 
In contrast to such in-plane inhomogeneity, the decoher- 
ence effect is small in SL's. Further, various conditions 
of the SL can be controlled by present experimental ca- 
pabilities. In fact, enhanced superconductivity in such 
systems was reported recently [HI, [3, [T3] . 

Here, we provide further insight into this problem by 
studying Hubbard and t-J models of SL's made up of UD 
and OD cuprate layers using layer extensions of the cel- 
lular dynamical mean-field theory (CDMFTjJfllll and 



the dynamic cluster approximation (DCA) [20|,l21| as well 
as the slave-boson mean-field (SBMF) theory. We find 
that the superconducting (SC) order parameter is indeed 
increased in UD/OD SL's as compared to its value in the 



uniform system. For some combinations of dopings, the 
SC order parameter is found to become larger than its 
largest value in the uniform system. Furthermore, the 
SBMF calculation predicts a transition temperature in 
the SL that exceeds the maximum T c in the uniform sys- 
tem. These results suggest that "higher- T c " supercon- 
ductivity may indeed be realized in SL's of cuprates. 

We first consider a quasi-two-dimensional (Q2D) Hub- 
bard model for the Cu d electrons: 

H H ub = E £{zi)n ia + n A n ii ~ ( c i<? c 3° + Hm ) 

i.o i (ij)£xy,a 

-f,£(4c,v + tf.c.). (1) 

Here, Ci a is the annihilation operator of an electron with 
spin a at site i and rii a — cl a Ci a , t(t z ) the nearest neigh- 
bor transfer integral along the xy{z) directions, U the 
on-site repulsive Coulomb interaction, and e(Zi) the on- 
site potential dependent on the layer index z%. In general, 
the potential s(zi) should be computed self-consistently 
by including the long-ranged Coulomb interaction be- 
tween electrons and between electrons and the back- 
ground charge. Instead of performing such a fully self- 
consistent calculation, we fix the potential profile as well 
as the Fermi level by focusing on short-period SL's which 
have a small number of inequivalent layers. Thus, the 
charge density profile should be regarded as a result of 
the charge transfer among the constituent layers @, 0] . 

To begin, we study the zero temperature SC order pa- 
rameters of SL and uniform systems using the CDMFT. 
We treat the interlayer correlations on the mean-field 
level [13] and take care of short-ranged intralayer cor- 
relations by using a square 2x2 cluster. In this approx- 
imation, the self-energy depends on the layer index z% 
and the intralayer cluster coordinates Rt. The CDMFT 
then maps the bulk lattice problem onto a number N z of 
independent effective cluster problems, with N z the num- 
ber of layers in the unit cell. The cluster problems are 
solved by using the Lanczos exact diagonalization (ED) 
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FIG. 1: (Color online) The superconducting order parameter 
$ as a function of carrier density n for the Q2D Hubbard 
model computed by CDMFT with the ED impurity solver at 
T = 0. Black symbols are the results for the uniform systems. 
Colored (shaded) symbols connected by a thin line denote the 
results for the SL. The inset shows the magnification around 
the optimally doped region. 



technique [HHi]. We use t/ = Wt and t z = 0.5*. Re- 
sults do not depend on the choice of these parameters in 
a significant way unless [/<ff (the total band width) 
or t z = (purely two-dimensional case) or t z — t (three 
dimensional case). 

Results for the <i-wave SC order parameter ^ = 
J2cr( ac iv c i+xcr) = - J2a {vc ia c i+ y S ) for the uniform sys- 
tem are shown as black dots in Fig. [1] A "dome" like 
shape is consistent with previous calculations for the two- 
dimensional Hubbard model [24J. We find an optimal 
carrier density n opt = 0.91 which maximizes VP. 

Next, we construct a SL in which UD layers with 
n > n opt and OD layers with n < n opt are repeated along 
the z direction. We plot ^f(zi) at the corresponding den- 
sities n = n(zi) in Fig. Q] as colored (shaded) squares 
connected by thin lines. Filled (open) squares are the 
results for SL's in which two (one) UD layers and two 
(one) OD layers are repeated, i.e. AABB (ABAB) super- 
lattices. One sees that VP in the OD layers is enhanced. 
The enhancement is stronger in AABB superlattices than 
in ABAB superlattces possibly due to the dilution in the 
latter. 

We also considered thicker SL's consisting of several 
UD layers and OD layers sandwiching nearly optimal- 
doped layers, e.g. AA'BCC superlattices. Even in this 
case, we found the enhancement in the OD layers. How- 
ever, in the nearly optimal-doped layers is suppressed 
because the pairing interaction is reduced by the OD lay- 
ers, and we never observed significantly larger order pa- 
rameters than in the uniform system. An example is 
shown by open circles in Fig. [1] 

To gain more insight into the mechanism responsible 
for the enhanced superconductivity in the SL and obtain 
estimates for T c , we next apply a SBMF approximation 



to a layer t-J model, the large U limit of the Hubbard 
model in which no double occupancy is allowed and near- 
est neighbor spins are coupled by exchange interactions 
J and J z . The Hamiltonian is given by 

H t j = ^e(zi)ri io - - f ^ (cj CT c> + H.c.) + J Sj ■ Sj 

i,a {ij)£xy,a {ij)exy 

-t, J2 (4*> + H - c ) + So , (2) 



where, Si = J2acr> Hcr^i^' with r the Pauli matrices, dia- 
ls a projected fermion operator defined as Ci a — Ci a {\ — 
nig-), and flic = c itJ Ci a . Here, we take t z = 0.5i, J = OAt 
and J z = for simplicity. The qualitative behavior does 
not depend on the choice of exchange interactions. 

In the SBMF approximation one introduces fcrmionic 
spinon f a and bosonic holon b operators with a lo- 
cal constraint b\bi + j\ a Ua = 1 an d rewrites 
the electron operators as Ci a = b\fi a . By in- 
troducing a mean-field decoupling in H t j; = 
Y,Mi*fi+*<r) = Y,Mi*fi+v°) [ umfo ™ resonance 
valence bond (RVB)], X { = EM&i+**)> A = 
E CT ( cr /^/»+^ff} = -E CT ( cr /^^+yff) (singlet RVB), X b = 
(b\b i+x ) = (b\b i+ y), xl = (Kh+z), and relaxing the lo- 
cal constraint to a global one, we obtain the mean-field 
Hamiltonians for spinons and holons 25, 26|, 27 1. 

Solving the self-consistent equations is straightforward. 
However, we find that the onset temperatures of x\ (and 
therefore xl) an d holon condensation Ao = (bi) 2 almost 
coincide for a relatively wide range of t z /t. Therefore, 
such a treatment is not useful to address the effect of 
the interlayer coupling on the SC transition. In fact, 
these phase transitions are artifacts of the mean-field de- 
coupling. To avoid these fictitious phase transitions, we 
replace x^ an d xl (weakly dependent on n and T at 
n <~ 1) in the holon Hamiltonian with a small number r 
of the order of 0.1, and perform a self-consistent calcula- 
tion for the other part. This is similar to approximating 
X b = S = 1 — n 25[. But, the present treatment can 
be used in a wider parameter range since the n and T 
dependences of x b an d x\ are included. There remain 
phase transitions associated with A and No. A higher 
transition temperature between them is expected to be- 
come a smooth crossover due to the gauge fluctuations, 
and only the lower temperature, i.e. SC critical temper- 
ature T c , remains as a true phase transition [281 ]. Since 
our main interest is in obtaining an estimate of T c , we 
keep our discussion on the mean-field level. 

Numerical result for the phase diagram is shown in 
Fig. [2] (a) as a function of carrier concentration n and 
temperature T. Below Tb, holon condensation occurs, 
and below T<\, A becomes finite. Thus, the lower tem- 
perature signals the onset of d-wave superconductivity. 
Optimal doping is located at n opt = 0.846 with a maxi- 
mum critical temperature T° pt = 0.132i. The character- 
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FIG. 2: (Color online) (a) Mean-field phase diagram of the 
Q2D t-J model as a function of carrier concentration n and 
temperature T with parameters as indicated. Below T, holon 
bose condensation No occurs, while below Ta the singlet RVB 
order parameter A becomes finite, (b) Order parameters as a 
function of n at T = O.Olt. Broken line: A, dash-dotted line: 
iVo, and solid line: d-wave superconducting order parameter 
9mf- The squares indicate the order parameters $mf of a 
SL in which layers with n = 0.7 and 0.9 alternate. 



istic T dependence of the order parameters of the uni- 
form system at under doping (n = 0.9) and over doping 
(n = 0.7) are displayed in the top panel and the bot- 
tom panel of Fig. [31 respectively. In Fig. [5] (b), we show 
the order parameters at low temperature as a function of 
n. In the UD region strong singlet correlations exist but 
weak coherence, while in the OD region this behavior is 
reversed. Consequently, the mean-field SC order param- 
eter >3>mf = N A has a maximum (*™/^ = 0.0451) as 
T c but at a slightly smaller carrier density n = 0.76. 

Next, we consider a SL in which UD (n = 0.9) and OD 
(n = 0.7) layers alternate along the z direction (ABAB 
stacking for simplicity). At low enough temperature, 
most holons are condensed giving rise to No ~ 0.1(0.3) 
in the UD (OD) layers. Since there are enough carriers 
(holons) to be paired in the OD layers, the proximity of 
A from the UD layers enhances the SC order parameter 
in the OD layer {^> MF = 0.0467 > as shown by 

the squares in Fig. [2] (b). On the other hand, in the UD 
layers, most carriers are already paired by strong pairing 
interactions. Thus, the interlayer coupling does not affect 
the UD layers significantly. These results are consistent 
with the CDMFT results for Hubbard SL's in Fig. [IJ 

The temperature dependence of the SL order param- 
eters is shown in the middle two panels of Fig. [3J The 
onset of A is dominated by the UD layer which has a 
higher Ta in the uniform system, so Ta = 0.14i. The 
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FIG. 3: (Color online) Mean-field order parameters of the 
Q2D t-J model as a function of T with the same parameters as 
in Fig. [2l Top (bottom): Uniform system with n = 0.9(0.7). 
Middle: Layers with n = 0.9 and 0.7 in the SL. Broken lines: 
A, dash-dotted lines: JVo, and solid lines: $mf(x5). Vertical 
line indicates T° pt — 0.132t of the uniform system. 



onset of No is roughly determined by the average carrier 
density n av = 0.8 and T, ~ 0.152i because it represents 
three-dimensional coherence. Thus, the lower tempera- 
ture corresponds to T c = Ta = 0.14i, which exceeds the 
maximum T c in the uniform system, T° pt = 0.132£, by 
about 6%. When the carrier density in the OD layer 
is further reduced, T, is increased but Ta is decreased. 
Therefore, an optimal combination is naturally expected. 

Finally, we check the finite temperature SBMF results 
using a layer extension of the DCA. The DCA for the 
SL also maps the bulk system onto a number of inde- 
pendent effective cluster models embedded in a dynamic 
host. In contrast to the CDMFT, the cluster is defined 
in reciprocal space. We solve the cluster problems with a 
noncrossing approximation (NCA) which has been shown 
to provide qualitatively similar results to a numerically 

~ 11. 



exact quantum Monte Carlo method [21 



As can be seen from Fig. [31 the order parameter in the 
SL system is very small close to T c . The same behavior 
is observed in the DCA-NCA calculations for the SL, 
which renders a reliable extraction of T c difficult if not 
impossible. We therefore present DCA-NCA results for 
the order parameter at a temperature below the T c of the 
uniform and SL system and argue that an enhancement of 
the order parameter in the SL close to T c should translate 
to an increase in T c over the uniform system. 

Figure [4] summarizes the results of the layer DCA cal- 
culations. Here we used t z — 0.5i, J = 0.3i, and an 
additional interlayer exchange J z — 0.25J. In the uni- 
form system we find an optimal doping n opt « 0.82 with 
T°P f w 0.0431 When rescaled by J, this T°p* is about a 
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FIG. 4: (Color online) The superconducting order parameter 
$ as a function of carrier density n for the Q2D t-J model 
computed by the DCA-NCA at a temperature T = 0.041t 
close to T c . The results for the uniform systems are indicated 
by the black symbols, while the colored (shaded) symbols con- 
nected by a thin line denote the result for the SL. Also shown 
is doping dependence of T c in the uniform system (open cir- 
cles). 



ers resulting from the proximity of strong pairing corre- 
lations from the underdoped layers. At certain combi- 
nations of doping, the superconducting order parameter 
even exceeds the largest value in the uniform systems. 
The slave-boson mean-field study shows a higher super- 
conducting transition temperature than the maximum T c 
in the uniform system. 
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factor 2.4 lower than that in the SBMF calculation. The 
order parameter 'F at a temperature T = 0.041i slightly 
below T c shows a similar n dependence as T c . For the 
AABB SL system, we again find a clear enhancement of 
the order parameter 'F in the OD layers which can even 
exceed the optimal value in the uniform system. This be- 
havior is consistent with the SBMF calculations at finite 
T and the CDMFT at T = 0. 

Calculations for t z < 0.5t (not shown) indicate that 
both the order parameter and T c are reduced, and, as a 
result, the enhancement of superconductivity. Note that 
in the limit of t z = 0, the enhancement has to vanish. 

The present mean-field calculations suggest a "recipe" 
to increase T c in SL's of cuprates: (1) Combine UD 
cuprates with large pairing scale and OD cuprates with 
strong coherence and (2) make the average carrier den- 
sity n av slightly smaller than the bulk optimal value n° pt . 
The latter condition is required to obtain a larger or- 
der parameter, which gives rise to higher T c even in the 
double-layer systems when there is enough phase stiff- 
ness. Requiring strong coherence, i.e., smaller n av than 
n opt , is related to the proposal in Ref. [irj(. A definite 
proof that this mechanism can give rise to higher T c re- 
quires calculations with fully three-dimensional clusters. 
Results of these calculations will be published elsewhere. 

To summarize, we studied the properties of model su- 
perlattices consisting of underdoped and overdoped high- 
T c cuprates. Two cluster dynamical-mean-field methods 
were extended for the superlattice geometry and applied 
to Hubbard and t-J models for the cuprate supercon- 
ducting supcrlattices and compared to a more transpar- 
ent slave-boson mean-field calculation. All calculations 
show enhanced superconductivity in the overdoped lay- 
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